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ABSTRACT 

High redshift galaxies selected on the basis of strong Ly-a emission tend to 
be young and have small physical sizes. We show this by analyzing the spectral 
energy distribution of nine Ly-a-emitting galaxies (LAEs) at 4.0 < z < 5.7 in 
the Hubble Ultra Deep Field . Rest-frame UV to optical 700A < A < 7500A 
luminosities, or upper limits, are used to constrain old stellar populations. We 
derive best fit, as well as maximally massive and maximally old, properties of all 
nine objects. We show that these faint and distant objects are all very young, 
most likely only a few million years old, and not massive, the mass in stars 
being « 10 6 — 10 8 M . Deep Spitzer Infrared Array Camera observations of 
these objects, even in cases where the object was not detected, proved crucial in 
constraining the masses of these objects. The space density of these objects, ~ 
1.25 x 10 -4 Mpc -3 , is comparable to previously reported space densities of LAEs 
at moderate-to-high redshifts. These Ly-a galaxies show modest star formation 
rates of ~ 8 M Q yr , which is nevertheless strong enough to have allowed them 
to assemble their stellar mass in less than a few million years. These sources 
appear to have small physical sizes, usually smaller than 1 kpc, and are also 
rather concentrated. They are likely to be some of the least massive and youngest 
high redshift galaxies observed to date. 

Subject headings: galaxies: evolution, galaxies: high redshift, galaxies: forma- 
tion, galaxies: structure, surveys, cosmology 
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Introduction 



The availability of the Infrared Array Camera (IRAC) of the Spitzer Space Telescope has 
recently made it possible t o measure the rest-frame optical light from high-redshift galaxies, 



up to redshiftsof ~ 6 (e.g.. lYan et al.l 



20051 : ISchaerer fc Pelld l2005t iDow-Hvgelund et al 



2006 



Chary et al. 



2005 



2005 



lEgami et al. 



20051 : lEvles et al. 



Stark et al.l 120071 ). Indeed, in just 



the past few years large numbers of Lyman break galaxies (LBGs) and Ly-a-emitting galaxies 
(LAEs) have been studied, and IRAC data have enabled mass and age estimates for these 
objects to be derived. Historically, LBGs and LAEs are two distinct classes of objects that 
have been observed over a very wide range of redshifts, thanks in large part to the fact that 
they are relatively straight-forward to identify by virtue of the presence of either a significant 
break in their spectral light distribution (LBGs) or the presence of a strong Ly-a emission 
line ( LAEs). As the number of known z > 5 LBGs and LA Es increased over the past few 



years (IDawson et al.l 12004 : lEvles et al.l 120051 : lYan et al.l 120061 . ; references therein) , the mass 



estimates of these objects have remained somewhat uncertain. Nevertheless observations 
have shown that some high redshift LBGs tend to already be quite massive (> 10 10 M ) and 
to be a few hundred million years old. The relatively large masses of these objects at a time 
when the universe was still quite young has been taken as an indication that these objects 
have had a star formation history (SFH) that allowed for most of their stellar mass to have 
formed early on. 

Mass estimates for LAEs, on the other hand, have been much more difficult to assess. 
Either the LAEs are very young, as indicated by the strong Ly-a line emission or the faintness 
of the continuum, for some other reason (e.g., extinction) makes it difficult to derive mass 



estima tes using their spectral energy distribution (SED) as is done for LBGs. iGawiser et al. 



( 120061 ). for example, required the use of a stack of 18 field LAEs to infer an average mass of 
pa 5 x 10 8 M Q . This low estimated mass, and the low luminosity of these objects in the rest- 
frame optical, is exactly what makes estimating their masses somewhat more difficult than in 
the case of LBGs. While the latter are usually selected based on broad band observations and 
colors, hence ensuring that these objects have bright continuum luminosities, LAEs selected 
based on the presence of Ly-a emission can be significantly fainter in the observed bands. 
For the handful of LAEs for which de tections by IRAC have been obta ined, a large mass 
in older populations has been derived (ILai et all 120071 : IChary et all 120051 ). There is however 
a detection bias here: the only published mass/age estimates so far are the ones for which 
IRAC detections exist; that is, they would tend to have significant old stellar populations. 

In this paper, we investigate the pop ulation of high redshift LAEs in the Hubble Ultra 
Deep Field (HUDF: lBeckwith et al.ll2007l ). a field that is relatively small (11 arcmin 2 ) but has 
the advantage of having been observed down to very faint magnitudes, and in a wide range 



- 3- 



of wavelengths. Our LAE sample makes no apriori assumptions as to the specific redshift of 
these sources, or as to whether these objects have been detected in all of the available data. 
We therefore aimed at reaching down to potentially fainter objects, and over a wider range 
of redshift, than previously done. Section [2] describes the object selection and the available 
data. In Section [3j we summarize the broad-band luminosities of these objects, while in 
Sections HJ [5j and [6] we discuss the morphology and SED-fitting process and provide mass 
and age estimates of the HUDF LAEs, irrespective of whether these objects were detected 
in the Spitzer IRAC bands or not. 



2. Observations 



The sources described in this paper were selected from the GRAPES su rvey (GRism 



ACS Program for Extragalactic Science, PI: S. Malhotra; see description in iPirzkal et al. 



2004 ). a slitless grism spectroscopic survey using the Advanced Camera for Surveys (ACS) 
on the Hubble Space Telescope (HST). The GRAPES project obtained low resolution spec- 
tra of 1421 objects in the HUDF down to a continuum level of z§5o = 27 mag. The 
GRAPES grism observations proved able to detect emission lines with fluxes as low as 
~ 5 x 10~ 18 ergs s _1 cm -2 in objects as faint as z 85Q = 29 mag in the conti nuum. A 
complete list of GRAPES emission line objects was published in Ku et al.l (120071 ) and con- 
tains 115 objects with redshifts 0.1 < z < 5.7. GRAPES identified a total of 9 LAEs at 
z > 4.0 (listed in Tabled]). We show individual spectra of these nine sources in Figure [TJ 
In most cases, the spectra shown are a narrow extraction of the GRAPES program data. 
In two cases, where the galaxy is comparatively large, we show a wider extraction (5225 
and 6139). Finally, in another two cases (631 and 9040), the GRAPES data have relatively 
few uncontaminated roll angles, and we therefore instead plot spectra from the subsequent 
PEARS program (Probing Evolution and Reionization Spectroscopically; PI: S. Malhotra), 
which reobserved the HUDF in four new roll angles. The emission lines in the largest objects 
(6139, 9040, and especially 5225) are not shown to best advantage in these one- dimensional 
extractions, because the spatial extent of the object gives a wide line-spread function in the 
slitless data, which tends to blur t he line. Detailed two dimensional images of the observa- 
tions of object 5225 are available in lRhoads et al.l (120051 ). The GRAPES emission line object 
selection was perform ed automa t ically using an emission-line search code and methodology 
described in details in lXu et al.l (120071 ). 



Some of the objects shown in Table [T] have rather faint emission lines. All lines were 
nevertheless significant at the level of 3.1 (N/3) 1 / 2 cr, where N is the number of GRAPES 
position angles in which the line was detected, and where N > 2 was required for inclusion 
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in ?. For an emission line object to be classified as a Ly-a galaxy, it had to also have (1) no 
other emission lines in the G800L grism wavelength coverage^] and (2) broad band photome- 
try in the HUDFimages consistent with a Lyman break at the redshift corresponding to the 
Ly-a line identification. Given less accurate photometry, a Ly-a line plus Lyman brea k 
may be difficult to distinguish from an [Oil] 3727A line plus 4000A break (IStern et al.ll2000l ). 
Fortunately, even the faintest object in our sample is detected witha signal-to- noise (S/N) 
greater than 5 two filters, and the rest are detected with (S/N)> 15 in two or more filters. 
This deep photometry shows that every object has a break where the flux is attenuated by 
a factor of at least 15 (see Table [3]). The strengths of these breaks are several times greater 
than would be expected from the D(4000) break, and are much more consistent with the 
decrement exp ected from absorption by the intergalactic medium (IGM) for objects at z ~ 5 
( IMadaulll995l ). Also, although it was not part of the GRAPES LAE selection process, these 
objects typically have very blue observed z — J colors (corresponding to ~ 1400 — 1800A 
rest-frame near-UV colors), based on the Near-Infrared Came ra and Multi- Object Spectrom- 
eter (NICMOS) HUDF observations jThompson et al.ll2005h . This makes it highly unlikely 
that any large amount of extinction is present in these objects and so further supports the 
interpretation of their colors as a Lyman break at high redshift. Further spectroscopy could 
certainly help secure the redshifts and Ly-a line identification for these objects. For exam- 
ple, higher resolution (R^j 1100") rmilH hfiln bv vprifvincr that the linps show the asvmmfitrv 

characteristic of Ly-a (e.g. 



Rhoads et al. 


2003; 


Dawson et al. 


2004) 



coverage could confirm the absence of other lines, and infrared coverage could confirm the 
presence of [Oil] A3727A at the same redshift as Ly-a . However, all present evidence points 
to these objects being Ly-a sources at z ~ 5.5. 

The HUDF, which li es within the Great Observatories Origins Deep Survey (GOODS; 
Giavalisco et al.l 120041 ) field, has the advantage of having been observed in a wide range 
of wavelengths. In this paper, we take advantage of the availability of HST NICMOS, VLT 
ISAAC, and Spitzer IRAC observations of this field, as well as the list of nine spectroscop- 
ically identified GRAPES LAEs. At the redshifts considered in this paper (4.0 < z < 5.7), 
combining these observations together allows us to examine these objects at rest-frame wave- 
lengths of approximately 700A < A < 7500A, and down to unprecedented depths. The ACS 
observations, using the B 435 , V 60 6, 1775, and z 850 bands e 



(Beckwith et al. 


2007 


). NICMOS 


• infrared ( 


Thompson et al. 


2005) 



using both the J and H bands. While the resolution of the NICMOS observations (« 0.2" 
pixel -1 ) is lower than that of the ACS observations, it is nevertheless high enough to clearly 



5 The case of Ly-a plus CIV A1549 emission did not occur in the GRAPES data set. 
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separate LAEs from their neighbors. We additionally used the ISAAC (Infrared Spectrom- 
eter and Array Camera) Ks-band observations of the HUDF, but the depth of these data is 
relatively shallow (25.1 mag) and with a resolution limited by atmospheric seeing (« 0.5"). 
Lastly, the HUDF was also observed in the infrared using IRAC as part of the GOODS 
Legacy program (M. Dickinson et al. 2007, in preparation). While the latter observations 
have the relatively low resolution of 1.2" pixel -1 , they offer a unique opportunity to probe 
the rest-frame optical properties of z ~ 5 LAEs. The IRAC data reach down to a depth of 
nearly 26 mag in both IRAC channel 1 (3.6 //m) and channel 2 (4.5 jttm). 



3. Photometry of Ly-a emitters 

The luminosity of each source was measured separately in each available band. Postage 
stamp images of the GRAPES LAEs are shown in Figure [2j Flux estimates were obtained 
using simple aperture photometry so that we could set meaningful upper limits to the flux of 
undetected sources, while the local noise level was estimated using a small annulus surround- 
ing each source. How we specifically dealt with non detections is explained in more detail in 
Section [51 We used nominally sized apertures, each containing as much of the object flux as 
possible while minimizing the flux contributions from nearby objects, by additionally using 
an object catalog derived using the ACS i^andz^o data. The latter allowed us to verify 
that contamination from neighboring sources remained small. 

Different aperture sizes were used in different bands. With the ACS data , we used an aper- 



ture size set to 5 times the size of the object as measured by SExtractor (IBertin fc Arnouts 



19961 ) (in this paper, "aperture size" refers to the diameter of the aperture), and no aperture 
correction was required. We compared these measured fluxes with SExtractor measurements 
and found that the two agreed very well, to better than 1%. We used an aperture of 1.25" 
with the NICMOS data and applied aperture corrections of 1.06 and 1.09 in the J and 
H bands, respectively (NICMOS Support 2007, private communication). When measuring 
fluxes from the ground based ISAAC data, we used an aperture size that was twice as large 
as the size of the source, as measured in the ACS z-band image and after accounting for 
a seeing disk of 0.5". As discussed above, the IRAC data both have a lower resolution 
and suffer from a higher level of contamination from neighboring sources than the rest of 
our data. The IRAC photometric measurement process was therefore more involved than 
when using ACS, NICMOS, and ISAAC data. Since many of the GRAPES LAEs were only 
marginally resolved in the IRAC observations, and because contaminating flux from nearby 
objects would cause flux measurements to be systematically higher (leading to systematically 
higher masses, as we discuss further below) , we developed a technique to remove as much of 
the contaminating flux as possible. Our approach consisted of performing two-dimensional 
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Table 1. GRAPES LAEs and Their Morphological Measurements. 



UID 


R.A. 

(J2000) 


Dec. 
(J2000) 


Redshift 

z 


Flux 

10~ 18 ergs s -1 cm~ 2 


M 

1500A 


-R50 
(kpc) 


C 


A 


631 


3:32:40.09 


-27:49:01.19 


4.00 


58. (20) 


-20.21 


0.67 


2.74 


0.097 


712 


3:32:42.81 


-27:48:58.51 


5.20 


17. (6) 


-21.31 


0.52 


2.53 


0.11 


4442 


3:32:39.38 


-27:47:39.17 


5.76 


24. (7) 


-19.54 


0.61 


2.89 


0.22 


5183 


3:32:34.51 


-27:47:27.11 


4.78 


32. (11) 


-19.82 


0.54 


2.84 


< 0.1 


5225 


3:32:33.26 


-27:47:24.76 


5.42 


22. a 


-22.77 


1.86 


2.18 


0.31 


6139 


3:32:37.95 


-27:47:10.99 


4.88 


60. (22) 


-22.14 


1.35 


3.19 


0.41 


9040 


3:32:41.08 


-27:46:42.45 


4.90 


36. (13) 


-22.41 


1.47 


2.65 


0.12 


9340 


3:32:40.67 


-27:45:56.11 


4.71 


49. (15) 


-20.52 


0.66 


2.69 


0.89 


9487 


3:32:40.17 


-27:46:00.55 


4.10 


55. (17) 


-21.98 


0.62 


2.36 


0.11 



Note. - The half-light radius (-R50), Concentration (C), Asymmetry (A) and absolute magnitude 
(Af^ oo O were estimated at a rest-f rame wavelength of 1500A. The redshifts and the line fluxes are from 



Xu et al.l ( 20071 ): the line fluxes from lXu et al.l ( 20071 ) shown in parenthesis were corrected using an aperture 
correction that we derived taking into account the size of each individual source. Units of right ascension 
are hours, miniutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 



tehoads et all |2005l ) 



- 7- 



general Sersic profile fits of all neighboring sources (whose p ositions were deri ved using the 
ACS i 775 andz 850 object catalog) using the program GALFIT (jPeng et al.l 120021 ). The center- 
most 4" region surrounding each LAE was masked and therefore not used during the fitting 
process. Each object was individually and interactively processed until the best possible fit 
was achieved, while avoiding over-subtraction of the flux contribution from the neighboring 
sources. Following this, the two dimensional models of the light profiles of nearby sources 
were subtracted from the original IRAC images and aperture photometry was performed on 
the resultant images, just as we did with the ACS, NICMOS, and ISAAC data. 

The reduced \ 2 °f the GALFIT fitting of the nearby sources was usually near unity. In 
the few cases where nearby sources were not perfectly fitted using a Sersic profile, we made 
sure that we did not end up over-subracting the amount of light contributed by these nearby 
objects. At the presumed redshifts of these objects, undersubtracting the contaminating light 
produced by nearby object ultimately results in an overestimate of the rest-frame optical 
flux of the GRAPES LAEs. This, in turn, can only lead to an inflation of the stellar mass 
estimates. We used a fixed aperture of 3" for all IRAC data and applied aperture corrections 
of 1.78 and 1.77 to IRAC channel 1 and channel 2, respectively (B. Mobasher 2006, private 
communication) . T he IRAC aperture cor rections were determined following extensive Monte 
Carlo simulations (IMobasher et al.l 120051 ). Note that we did not use IRAC channels 3 and 
4, as these proved a bit too noisy for our purpose. Table [3] shows the magnitude estimates 
for each source and in each available band, as well as the 3 a upper limit estimates for non 
detection. 



4. Morphology 



We measured the size and morphology of each GRAPES L AE using both GALFIT 
( IPeng et al.l 120021 ) and the CAS parameters (jConselice et al.l 120001 ) . We computed the rest- 
frame size, concentration (C) and asymmetry (A) of each source separately in each of the 
available ACS bands. 1500 A rest-frame values of these parameters were then computed by 
interpolating these measurements. 

The individual measurements are shown in Table [TJ GRAPES LAEs are found to be sharply 
concentrated and asymmetric sources, with average values of 2.67 ± 0.30 and 0.23 ± 0.30 
for C and A, respectively. The half-light radii obtained using GALFIT are consistent with 
this, with an average size of 0.92 ± 0.50 kpc (or 0.17 ± 0.07") while best-fitting light profiles 
were sharp and concentrated in most cases. Three of the sources (5225, 6139, and 9040) are 
significantly larger than the rest of the GRAPES LAEs, however, with an average half light 
radius of 1.56 ± 0.27 kpc, as shown in Figure[3l Moreover, two of these three objects (5225 



Table 2. GRAPES LAEs Photometric Measurements 



UID 




ACS 




NICMOS 


ISAAC 


IRAC 


B435 


V606 


1775 


z 850 


J1100 


Hl600 


Ks 

2200 A 


IRAC 

3600A 


IRAC 

4000A 


631 


29.59 


OR 87 


OR A 1 
ZO.41 


26.47 


OR AO 

zo.4y 


0(\ 7fi 
ZO. / D 


> 26.61 


07 1 7 
Z f . 1 / 


97 1 R 
Z / .10 




(0.35) 


/a no\ 

(0.02) 


/A A 1 \ 

(0.01) 


(0.03) 


(0.07) 


(0.10) 


(0.83) 


/A 0^7\ 

(0.37) 


/A A A \ 

(0.44) 


712 


31.48 


90 RR 


97 9^ 
Z / .ZD 


27.12 






26.74 


9& 

Zo. uo 


^ 97 A 
> Z t .40 




(2.01) 


(0.25) 


(0.03) 


(0.04) 






(1.32) 


(0.84) 


(0.83) 


4442 


> 32.96 


31 RA 
01. 04 


90 

zy .oo 


28.36 


90 /1R 


0U.O0 


> 27.78 




-v. 90 n7 




(0.83) 


(1.22) 


(0.18) 


(0.11) 


(0.52) 


(2.04) 


(0.83) 


(0.83) 


/A O O \ 

(0.83) 


5183 


31.76 


28.27 


27.39 


27.97 


28.60 


28.11 


> 27.33 


> 28.13 


> 27.85 




(2.04) 


(0.06) 


(0.03) 


(0.09) 


(0.28) 


(0.20) 


(0.83) 


(0.83) 


(0.83) 


5225 


28.88 


28.29 


25.97 


25.93 


26.05 


26.68 


26.54 


26.51 


25.86 




(0.41) 


(0.17) 


(0.02) 


(0.04) 


(0.07) 


(0.16) 


(1.62) 


(0.28) 


(0.26) 


6139 


> 32.92 


26.93 


25.52 


25.66 


25.83 


25.84 


26.20 


26.35 


27.10 




(0.83) 


(0.04) 


(0.01) 


(0.02) 


(0.12) 


(0.15) 


(0.90) 


(0.23) 


(0.39) 


9040 


30.14 


28.30 


25.82 


26.23 


26.16 


26.11 


27.06 


> 28.45 


> 28.36 




(1.29) 


(0.17) 


(0.02) 


(0.05) 


(0.14) 


(0.22) 


(2.49) 


(0.83) 


(0.83) 


9340 


33.04 


28.09 


27.51 


27.69 






> 27.40 


> 27.48 


> 27.16 




(0.83) 


(0.05) 


(0.03) 


(0.07) 






(0.83) 


(0.83) 


(0.83) 


9487 


30.07 


27.25 


27.16 


27.27 






> 27.24 


27.47 


> 27.93 




(0.47) 


(0.02) 


(0.03) 


(0.05) 






(0.83) 


(0.47) 


(0.83) 



Note. - 



Error estimates are shown in parentheses; 3 a upper limits are indicated with greater-than signs. 
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and 6139) have larger measured asymmetry values than the rest of the sample. In general, 
objects with larger asymmetries are observed to be have large half-light radii, raising the 
possibility that these objects could in fact be the result of smaller interacting components. 
The only exception is object 9340 which, while small, is also quite irregular in shape. The 
estimated rest-frame 1500 A siz es are consistent with the expected siz es of objects at thes e 
redshifts as shown in Figure 1 of iFerguson et al.l (120041 ) and Figure 7 of iPirzkal et al.l (120061 ). 
As expected if the size of these objects scale approximately as H" 1 , the GRAPES LAEs are 
indeed smaller than what was observed at lower redshifts by [Ferguson et al.l (12004 ). On the 
other hand, we found that the mean ellipticity an d concentration o f the GRAPES LAEs are 
both marginally lower than what was measured bt IFerguson et al.l (120041 ) and shown in their 
Figure 3. The mean GRAPES LAE 1500A ellipticity and concentration are observed to be 
0.31 ± 0.18 and 2.67 ± 0.30, respectively. 

Morphologically speaking, we observed GRAPES LAEs to be small, concentrated sources 
with somewhat irregular shapes, reminding us of local dwarf galaxies. The small sizes of the 
GRAP ES LAEs are also similar to the sizes reported for z ~ 6 LBGs by lDow-Hygelund et al. 



( 120071 ) . Several sources (5225, 6139, 9040, 9487) show a resolved clumpy or complex structure 
(44% of all sources), while six objects have an effective radius -R 50 < 1. 3kpc (67% of all 



sources). For purposes of comparison with sources at a lower redshift, IVenemans et al. 
( 120051 ) found that, based on 40 z ~ 3 sources, 19% of their sample was observed to have 
clumpy structure, and 25% of their sources had an effective radius R^o < 1.3kpc. Even 
considering the relatively low number involved, GRAPES LAEs appear to be systematically 
smaller, and maybe more irregularly shaped than LAEs at the lower redshift of z — 3.1. 



5. Population Synthesis 



The magnitudes shown in Table [3] provide us with unprecedented measurements of 
the rest-frame UV-to-optical light distribution of these faint high-redshift objects. A first 
glance, Table [3] shows these objects to be brighter in the i 775 , z 850 , Jnoo, and H 1600 bands 
(corresponding to rest-frame 1200 A < A < 2700 A) than in the IRAC bands (optical rest- 
frame wavelength of 6000A < A < 7500A) and than in the B 435 and V606 bands (rest-frame 
UV of 700 A < A < 11 00 A). The lack of strong emission in the B 435 and V606 bands is due 
to HI absorption from the IGM of these objects blue ward of the Ly-a limit. 

It is worth emphasizing again that all LAEs were selected spectroscopically and that 
these sources were not selected using broad band colors. They nevertheless all exhibit very 
blue intrinsic colors, as shown in Fi gure HI and app ear t o be signifi c antly bluer than the 
LAE samples previously observed by lLai et al.l (120071 ) and lEyles et al.l (120051 ). One is imme- 
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diately drawn to the possibility that these objects might well contain a population of active, 
young stars, which would of course be consistent with the detection of Ly-a emission in 
their GRAPES spectra. There are no known X-ray emitters wi thin 5" of any of the sources 
presented here in the Chandra lMs observations of the field (jGiacconi et al.ll2002l ). which 
excludes the blue nature of these objects from being caused by an active galactic nucleus, 
confirming earlier conclusions of Malhotra et al. (2003), Wang et al. (2004), Gawiser et al. 
(2006). 

A more quantitative determination of the age of these objects, or at least of the age of the 
stellar populations of these objects, as well as an estimate of their stellar mass, can be ob- 
tained by comparing the spectral energy distribution of these sources with stellar population 
synthesis models with known star formation history, masses, metallicities, and extinction 
values. 



We compared the GRAPES LAEs' SEDs with models from lBruzual fc Charlotl (120031 . here- 
after BC03), using the Paduva evolutionary tracks and using a Salpeterl ( 1955 ) initial mass 
function (IMF). The choice of a Salpeter IMF allows us to directly compare GRAPES LAE 
mass estimates with ones previously published in the literature. It also allows us to view 
our results as upper limits to the stellar mass of the GRAPES LAEs: because the Salpeter 
IMF continues as an unbroken power law to low stellar masses, it can overestimate the stel- 
lar masses of objects by a factor of 2 or more if there is a turnover at ~ 1 M (see, e.g., 



Baldrv fc Glazebrookl 120031 : Bell et al.ll2003l : lLai et all 120071 ). The initial goal of this effort, 



and since many of the magnitudes listed in Table [3] are either upper limits or have large 
associated errors, was to determine whether some reasonable constraints could be set on 
some of the physical properties of these objects; that is, we might not be in a position to 
unambiguously determine their exact nature, but we might be able to determine what they 
are likely not to be. 

We used three different set of BC03 models, each with a different SFH: (1) a single instan- 
taneous burst of star formation (SSP), (2) a single but exponentially decaying burst, with 
an e-folding time of r (EXP), and (3) two separate instantaneous star formation bursts sep- 
arated by an arbitrary amount of time, tb U rst ; an d with each burst contributing an arbitrary 
fraction of the final stellar mass of the galaxy (2BP). When using either of these three SFHs, 
the only free parameter that was allowed to vary during the fitting process was the total 
stellar mass of the model. The redshift, metallicity, age, and extinction were not treated as 
free parameters but as input assumption of the models. 

Using each of these three SFH scenarios, we generated an extensive set of simulated spec- 
tra. Model spectra were computed at 27 discrete ages, ranging from as low as 5 Myr and 
up to the age of the Universe at the redshift of interest (^ 1.2 Gyr; we adopt a ACDM 
cosmology throughout this paper with Qm = 0.3, Q\ = 0.7, and H = 70). The redshifts 
were fixed to the spectroscopically determined redshifts listed in Table [TJ For the SSP and 
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EXP SFH scenarios, we generated model spectra using several values of the metallicity Z 
(0.0001, 0.0004, 0. 004, 0.008, 0.02 [ZJ , and 0.05) and extinction values (0 < A v < 3, using 
the Calzetti law, ICalzetti et al.ll2000i ). The 2BP models assumed a more restricted range 
of metallicities, using a low metallicity for the first burst (Z=0.0001) and a solar metallicity 
(Z=0.02) for the second burst, but allowed for the same range of extinction values as the 
SPP and EXP scenarios. Differential extinction was not included in the 2BP model. We 
recognize the fact that the two stellar populations could in principle be subject to different 
amounts of dust extinction. However, GRAPES LAEs are observed to be intrinsically very 
blue objects and hence to be objects that are dominated by a combination of young stars 
and low extinction. The data are not sufficient to allow us to study in details whether these 
two stellar popul ations do in fact suffer from slightly different amount of dust extinction. 
IGM absorption (lMadaulll995l ) was added to all of the model spectra. A total of ~ 80,000 
separate model spectra were generated at each of the nine redshifts shown in Table [T] 



Considering each GRAPES LAE separately, we began the fitting process by minimizing 
the error-weighted x 2 P er degree of freedom (xt) between observations and a particular 
model spectrum. During this step, a model spectrum was only allowed to be scaled up or 
down by an arbitrary amount, thereby determining the mass of the galaxy assuming that 
the model spectrum was a good fit to the observations. We only used bands where an object 
was detected, and in addition, we excluded observations in the 1775 band, which we knew to 
be contaminated by potentially a large, but uncertain, amount of Ly-a emission. We did 
not attempt to simply correct the broad band 1775 measurements using the fluxes listed in 
Table [lj because of the large uncertainties associated with the aperture correction applied 
to the measured slitless spectroscopic fluxes. We also excluded the B 435 and V606 bands that 
we knew to be strongly affected by IGM absorption. 

This fitting process, while determining a mass estimate for each object, does not however take 
into account any of the upper limit estimates in cases of a non-detection. We incorporated 
this additional information in two dist inct ways. First , we computed the Kendall rank 
correlation value for censored data sets (llsobe et al.lll986l ) between individual models (each 
scaled to best fit bands where the object was detected) and the observations. This method 
allowed us to quantitatively determine how badly upper limit constraints might be violated 
by a particular model, but it has the drawback of not being able to take into account the 
errors in the flux measurements. It has the additional drawback of producing a goodness-of-fit 
estimator that is strongly quantized because of the limited number of available measurements. 
While models that strongly violated upper limit constraints could be identified, many of our 
models led to similar Kendall correlation values, making it difficult to select one model over 
another. We therefore also relied on another method based on a more generalized version of 
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the error-weighted xl described above. This method takes into account upper limit estimates 
and assigns a xl penalty when a model predicts flux values that exceeds established upper 
limit values (xl c )- The value of this additive penalty is — lnl, where I = P(< cr|/x) is the 
probability that the actual flux is less than a if the flux is assumed to be normally distributed 
with mean \i and width a, where a is the flux detectio n limit of t he da ta and \i is the flux of 



the model in that bandpass. This method was used in lLai et al.l (120071 ) where it is described 
in more details. 

We found a good agreement between the Kendall rank method and the xlc method, 
and both methods could indeed be used to reject models that over predicted fluxes in non 
detected bands. However, and because the xlc method allowed us to take flux measurements 
into account, which was crucial for some data points (such as the relatively low signal-to- 
noise ratio ISAAC Ks observations), we eventually relied more on the xl c method. Models 
were thus rejected by first computing individual xlc values and then subsequently determin- 
ing whether a given xlc value was an indication that a model was a statistically bad fit to the 
observations. The exact cut-off value of xlc use d to reject a set of models was determined 
separately for each object since the number of bands used to compute xlc differed from one 
object to the next. We defined P(xl c > xl co ) as ^ ne probability of obtaining a value of xl c as 
large as the observed xl co ^ the observations really did follow the assumed distribution (see 



Taylor! 1 199 71 ). We then assumed that observations and models differed significantly when 
Probd(xlc > xlco) was ^ ess than 5% (2a level) and rejected the corresponding model on the 
basis that it was a poor fit to the observations. Models that could fit the observations could 
often be identified, using one of the three SFH scenarios that we considered. Some notable 
exceptions is object 5183, that we found to be somewhat poorly fitted by all of the models 
that we generated. Cut-off values of P(xl c > xl Co ) were therefore loosened for these two 
objects, and hence the selected best fits for these objects have higher xlc values than the 
rest of the GRAPES LAEs. 



It is important to stress that we found that properly accounting for the non detection of 
an object, particularly in the IRAC bands, was crucial when attempting to either include or 
exclude a particular set of models. The IRAC fluxes, especially the ones that we corrected for 
the effect of neighbor contamination, allowed us to determine an upper limit to the optical 
rest frame flux emitted by the source. We found that such an upper limit imposed a strict 
upper limit to the number of older stars that could be accounted for in that object. This in 
turn allowed us to exclude a large subset of our original models. In the following section, 
we examine the mass estimates for the GRAPES LAEs and further summarize the range of 
acceptable input BC03 model parameters. 
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6. Mass and Age 

As expected with such a large number of input parameters, we found that a wide range 
of BC03 models provided an adequate fit to the observations. While a best fitting model 
could always be identified by simply selecting the model with the lowest value of xl c (which 
we call the best fit model), we also determined the largest mass estimate and the oldest 
age estimate for each GRAPES LAE by simply examining the subset of BC03 models, and 
hence the corresponding parameters that went into creating these models, that were not 
rejected at the 2a level. The model producing the best fit to observations was not always 
the most massive, or the oldest, model that properly fitted the observations. We defined 
the maximum-mass models and the maximum-age models as those with the largest mass 
and age, respectively, that were good fits to the observations (e.g., 2 a level or higher when 
possible, or highest possible confidence otherwise). 

Table [3] summarizes our results using the SSP SFH scenario. This is a simple model 
that assumes that all of the stars were created in a single burst event and then simply 
passively evolved. Examining the best-fit models shown in Table |3j we find that all LAEs 
appear to have very young stellar populations (a few times 10 6 yr) and low stellar masses, 
with only one object (6139) having a mass estimate that exceeds 10 8 M Q . Extinction values 
were consistently low, while the metallicities were found to be subsolar in most cases. How- 
ever, it should be noted that metallicity was not well constrained and that a wide range of 
metallicities were found to be acceptable. Naturally, a wider range of models were found to 
properly fit the observations and the masses and ages of the GRAPES LAEs could possibly 
be somewhat larger than those of the best-fit models, as shown on the last two columns of 
Table [31 Similarly, a range of age estimates were found to be acceptable. Yet, we found that 
acceptable ages remained below a few tens of millions years, while masses remained below 
a few times ~ 10 s M . Note that none of our BC03 models successfully fitted object 5183, 
ot 9040 [P(Xvc > Xvcq) < and that all of our BC03 models could hence formally be 
rejected at the 2 a level. In these cases, Table |3] shows the model with the highest possible 
confidence level. 

Table H]is similar to Table [3]but shows the results obtained using the EXP SFH scenario. 
In this scenario, star formation occurred in a burst that was not instantaneous but instead 
exponentially decayed as a function of time, with an e-folding time of r. Unlike the previous 
SSP model, not all of the stars have the same age in this scenario, and this scenario allows 
for an increasingly larger proportion of stars to be older as the galaxy ages. The EXP models 
would therefore be expected to result in slightly larger age and mass estimates than the SSP 
scenario. As shown in Table HJ this is indeed marginally the case. The maximum ages shown 
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in Table H] are found to be in some cases up to 10 times larger than those produced using 
the SSP scenario. The maximum masses, in most cases, did not exceed a few times 10 8 M 
and were only slightly larger than the maximum masses derived using the SSP model. 

The mass and ages estimates derived using our third SFH scenario (2BP) are shown 
in Table [5j By design, the two-burst scenario was conceived to allow us to fit the SEDs of 
objects that are intrinsically very blue in their rest-frame UV light yet appear to produce 
a significant amount of light in the rest-frame optical. The 2BP scenario allows for the 
presence of two distinct stellar populations, one young and blue and another older and 
therefore relatively redder. This scenario should therefore provide better fits to observations 
with both large ACS and IRAC fluxes. It should also be expected to yield higher mass 
esti mates in cases o f high upper flux limits in the IRAC bands. This scenario is similar to 



the iDickinson et al.l (120031 ) maximally old population model. As shown in Table the age 
estimates derived from the best fit to the data are in some cases significantly larger than 
those derived using either the SSP or EXP scenario. This is of course achieved by allowing 
a significant fraction of the stellar population of the galaxy to be in the form of a maximally 
old (> lGyr) population of stars. As previously observed, extinction values were low, and 
we found that lower extinction values were generally preferred using this SFH scenario. But 
we did find that a range of extinction values and burst mass ratios could fit the observations, 
leading to larger mass estimates, as long as at a few times 10 6 M Q of the galaxy mass was 
in the form of a young stellar population with an age of ~ (5 — 35) x 10 6 yr. The younger 
stellar population accounts for the significant UV flux in these objects, while the much older 
stars (often up to the age of the universe at the object 's redshift) account for an increase of 
the total mass of the galaxy by a factor of about 2, when compared with the maximum- mass 
SSP results. 

Figure shows the best-fit 2BP model for each GRAPES LAE (solid lines), as well as 
the maximum-mass selected from models that were not rejected at the 2 a level (dashed 
line). Note that, as mentioned above, two of the GRAPES LAE sources (5183, and 9040) 
are not very well fitted by any of our chosen models. It is worth addressing again the possi- 
bility that one, or both objects could be a low redshift interloper. However, and as pointed 
out in Section [21 every single one of these three objects has a very significant color break, 
consistent with what we would expect from IGM absorption at high redshifts, as well as 
very blue observed z 850 -Jnoo colors. We performed additional simulations of the B 4 3 5 -z 8 5o 
and z 85 o-Jiioo colors of these objects that would be expected if they were at much lower 
(z Rj 0.25 — 0.5) redshifts, assuming that the detected line was not Ly-a but [OIII] A5007A, 
and that the color break was the result of the 4000 A break. We concluded that the presence 
of a strong 4000 A break, together with a very blue rest-frame near-UV color, cannot be 
properly accounted for anything but a young stellar population at high redshift, indepen- 
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dently of the amount of extinction present, the metallicity of the stellar population, and the 
age of the stellar population. The latter led us to favor the interpretation that objects 5183 
and 9040 are indeed very faint high redshift Ly-a emitters. We recognize however that fu- 
ture high-resolution spectroscopy will be required in order to unequivocally address this issue. 



7. Discussion 



Based on the detection of nine Ly-a sources in the HUDF, and using a 1/V max method- 
ology, we derive a mean number density of ~ (1.25 ± 0.6) x 10 -4 Mpc -3 . Comparison 
with previous studies of high redshift Ly-a sources is difficult, since each study has slightly 
different limiting Ly-a line fluxes and luminosity selection biases. This work in particular 
does not claim to be a complete sample of LAEs at these redshift s. Still, down to a limiting 



Ly-a line luminosity of 2.0 x 10 42 ergs s 1 . iTaniguchi et al.l (120051 ) found that the density of 



z ~ 6.6 Ly -a sources is 1.2 x 10 -4 Mpc -3 . They furthermore estimated, based on observa- 
tions from lOuchi et al.l (120031 ) . that the density of z = 4.9 Ly-a source s is 2 x 10~ 4 Mpc -3 , 
down to the same line luminosity limit. At somewhat lower redshifts, iGawiser et al.l (120061 ) 
recently estimated that the number density (again down the same limiting line luminosity 



limit ) of z = 3 



Ly-a sources is (3 ± 1) x 10 4 Mpc 3 . If we restrict our sample to the 



Taniguchi et al.l (120051 ) limiting line luminosity, we still detect all 9 sources, and the HUDF 
LAE density remains ~ (1.25 ±0.6) x 10 -4 Mpc -3 , which is consistent with previous works. 



Simila rly, we can also compare our observations with the results from iRhoads fc Malhotra 
(120011 ) at z = 5.7. The latter observations were car ried out down to the brigh ter line lu- 
minosity limit of 5.3 x 10 42 ergs s -1 . In their work, IRhoads fc Malhotra! (120011 ) inferred a 
z = 5.7 Ly-a density of ~ 4 x 10 -5 Mpc -3 . Down to their brighter line luminosity limit, we 
detect seven LAEs (631, 4442, 5183, 6139, 9040, 9340, and 9487) and der ive a Ly-a density of 
(9.68±4.8) x 10 -5 Mpc -3 . Similarly, in this case ITaniguchi et al.l (120051 ) comput ed a. z ~ 6.6 
Ly-a object density of 2 x 10 -5 Mpc -3 , while the z ~ 4.9 object density from lOuchi et al. 
( 120031 ) becomes 4 x 10 -5 Mpc -3 . Our observations therefore co nfirm the relatively low spac e 



density of Ly-a objects at these high redshifts, as seen also by iMalhotra fc Rhoadd (120041 ) 
confirming an apparent lack of, or at least slow, redshift evolution of the density of LAEs. 



We found that, no matter which SFR scenario we considered, the GRAPES LAEs were 
unlikely to be very massive. Both the single instantaneous burst (SSP) and the exponen- 
tially decaying burst (EXP) scenarios resulted in similarly low mass estimates, the best EXP 
models appearing to be those with low values for the r-parameter. The double-burst see- 



Table 3. Best Single Instantaneous Burst SFH (SSP) Fits to the Observations. 



UID 




Mass 


Av 


z 


PCy 2 > v 2 1 




Mass 




10 6 yr 


10 8 M Q 




%z Q 




10 6 yr 


10 8 M Q 


631 


4.5 


1.27 


0.4 


2.0 


70.21 


45.0 


5.0 


712 


3.5 


0.3 


0.0 


100.0 


93.21 


45.0 


13.2 


4442 


2.5 


0.05 


0.0 


0.5 


96.95 


200.0 


5.14 


5183 a 


2.5 


0.90 


0.6 


2.0 


3.46 


2.5 


0.98 


5225 


4.5 


1.34 


0.0 


0.5 


56.30 


9.0 


2.89 


6139 


20. 


10.0 


0.4 


0.5 


78.91 


30.0 


29.0 


9040 a 


2.5 


0.89 


0.0 


0.5 


0.05 


2.5 


0.94 


9340 


0.5 


0.16 


0.0 


250.0 


100.00 


20.0 


2.57 


9487 


1.0 


0.40 


0.05 


20.0 


6.46 


9.0 


0.82 



Note. — Listed are the age, mass, extinction, and metallicity that resulted in the 
lowest values of xl c - 

a Object failed to be properly fitted at the 95% confidence level 



Table 4. Best Exponentially Decaying Star Burst SFH (EXP) Fits to the Observations 



UID 


A crp 


Mass 




Av 


z 


P(y 2 > v 2 1 


A cp 

A 6 c max 


Mass 




10 6 yr 


1O 8 M 


Gyr 




%z 




10 6 yr 


1O 8 M 


631 


8.0 


0.65 


0.005 





250.0 


71.06 


300.0 


11. 


712 


4.5 


0.30 


0.005 


0.0 


250.0 


94.09 


300.0 


14.8 


4442 


2.5 


0.07 


0.001 





0.5 


96.84 


1200.0 


6.76 


5183 a 


3.0 


0.81 


0.005 


0.6 


2.0 


3.38 


3.0 


0.81 


5225 


5.0 


1.41 


0.001 


0.05 


0.5 


54.41 


300.0 


25.3 


6139 


20. 


18. 


0.001 


0.4 


0.5 


80.54 


150.0 


27.0 


9040 a 


2.5 


1.2 


0.001 


0.0 


2.0 


0.03 


2.5 


1.24 


9340 


1.0 


0.34 


0.001 


0.0 


0.5 


100.00 


100.0 


2.63 


9487 


3.0 


0.32 


0.0125 


0.0 


20.0 


6.67 


15.0 


0.66 



Note. — Listed are the age, mass, extinction, and e-folding time r that resulted in the 
lowest values of xl c - 

a Object failed to be properly fitted at the 95% confidence level 



Table 5. Best Double Instantaneous Burst SFH (2BP) Fits to the Observations 



UID 




Mass 


2nd Burst 


Av 


Pf v 2 > Y 2 ) 
VAi/c ^ Ai>cq) 




Mass 




10 6 yr 


10 8 M Q 


Age (10 6 yr) 


Mass (%) 






10 6 yr 


10 8 M Q 


631 


8.5 


1.56 


8 


50 


0.15 


73.42 


1209.0 


30.9 


712 


16. 


0.68 


3.0 


50 


0.2 


94.31 


1010. 


22.7 


4442 


2.0 


0.065 


1 


1 





96.82 


870.0 


11.4 


5183 a 


17. 


1.09 


13. 


1 


0.2 


3.28 


17. 


1.08 


5225 


4.5 


1.28 


0.5 


2 


0.0 


56.42 


1000. 


73.8 


6139 


20. 


8.70 


4.5 


10 


0.2 


81.28 


1006. 


50.38 


9040 a 


7.0 


1.36 


2.5 


1 


0.0 


0.04 


7.0 


1.36 


9340 


1.5 


0.18 


1.0 


99 


0.0 


100.00 


1200. 


7.47 


9487 


1.0 


0.61 


0.5 


8 


0.15 


6.30 


1200. 


1.85 



Note. — Listed are the age, total mass, and extinction of the object, as well as the age of the second 
burst and the fraction of the total stellar mass it contains, that resulted in the lowest values of \l c - 

a Object failed to be properly fitted at the 95% confidence level 



-19- 



nario (2BP) led to mass estimates that were on average twice that of the other two SFH 
models. The latter was expected, since the 2BP scenario was conceived to allow fo larger 
stellar mass estimates by asking the question how large a fraction of the stellar popula- 
tion of a GRAPES LAE could in fact be in the form of an older stellar population, with a 
lower mass-to-light ratio. The different mass estimates derived using our three distinct SFH 
models should be kept in mind when comparing the results presented here with those in 
the literature. iFinkelstein et all (120071 ) who observed narrow-band-selected Ly-a sources as 
part of the Large Area Lyman Alpha (LALA) survey, derived masses for 22 of their sources 
and estimated them to be in the range 2 x 10 7 to 2 x 10 9 M . They used the best fit of 
an exponentially decaying burst model (similar to our EXP SFH model). The masses we 
derived for the GRAPES LAEs using th e best-fitting EXP SFH models (Table g) are about 
a factor of 2-3 lower than those found by IFinkelstein et all (120071 ). with masses varying from 
0.1 x 10 8 to 17. x 10 8 M . The largest mass-estimates shown in Table H] do not exceed 
2 x 10 9 M . Using recent observati ons of z « 5.7 Ly-a sources in the GOODS-N field and a 
simple SSP model, lLai et all (120071 ) estimated the masses of their t hree sources to be in the 
range (1.4 — 5) x 10 9 M Q , consistent overall with the results from IFinkelstein et all (120071 ) 
but, again, significantly larger than most of the best fit Ly-a mass estimates, as well as many 
of the maximum Ly-a mass estimates derived using the SSP, EXP, or even 2BP SFH model 
(Tables EH and EJ. 

These differences in mass estimate are unlikely to be caused by our respective choices of 
BC03 templates. Indeed, they are most likely the results of our reach to fainter luminosities 
(by nearly 2 AB magnitudes) and hence masses. Our sample was not selected based on the 
presence of a strong broad band continuum emission and hence might constitute a more 
complete sample (at least in term of masses) of LAEs. As they stand, the GRAPES LAEs 
are objects that appear to be significantly less massive than previ ously observed LA Es at 
high redshift, as well as LBGs such as the i-dropout population from Eyles et al. ( 20051 ) [with 
estimated masses of ~ (1.3 — 3.8) x 10 10 M Q ]. This is consistent with these objects being 
observed to produce most of their light in the rest-frame UV (i.e. a young stellar population 
with a relatively lower mass-to-light ratio) and to be either faint (e.g., objects 712 and 9487) 
or even undetected (e.g., objects 4442, 5183, 9040, and 9340) in the rest-frame optical, where 
one would expect light produced by an older, higher mass-to-light ratio stellar population. 



We can furthermore examine the star formation rate (SFR) of the HUDF LAEs using 
either their Ly-a line luminosities, listed in Table [TJ or their rest-frame UV luminosities 
(jKennicuttlll998l ). We can additionally use the SFR estimates that the BC03 code produces 
when generating models. Using the Ly-a fluxes listed in Table (U we derive SFRs that are 
consistently low for all sources, with an average SFR of 7.8 ± 3.2 M Q yr _1 . Similarly, the 
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BC03 SFR estimates imply an equally and systematically low SFR of 3.6±2.0 M Q yr _1 . SFR 
estimates derived using observed broad band rest-frame 2500A UV are much more uncertain 
and result in a higher average SFR of 18.9 ± 22.9 M Q yr _1 . However, and independently of 
which of these SFR estimates we consider, the low mass estimates we derived for the HUDF 
LAEs imply that they have had ample time to form their stars by the observed redshifts, 
even with an SFR that is only a few solar masses per year. Strictly speaking, there is no 
need for the SFR o f these objects to have been higher in the past, as was the case for the 
objects observed by lYan et al.l (120061 ). Higher SFR values would in fact cause these objects 
to have assembled their stellar masses even more quickly. Based on the observed SFR alone, 
the HUDF LAEs should have started to form stars no earlier than a few million years before 
we observed them, unless their SFR was significantly lower in the past, consistent with the 
fact that we have derived low ages for these objects. 

We are likely to have found very young and very low mass Ly-a emitters that formed recently 
as a result of reporting on a whole sample of Ly-a sources, including objects that were not 
detected in the rest-frame optical bands. The IRAC observations, when contamination from 
neighboring objects is properly accounted for, provide stringe nt limits to the maxim um stel- 
lar mass that these objects could have. It is noteworthy that iGawiser et al.l (120061 ) inferred, 
using a stack of 18 Ly-a sources at z ~ 3.1 Ly-a , that Ly-a sources at this redshift should 
have an average mass that is on the order of 5 x 10 8 M , using models with constant SFH 
(equivalent to our EXP models with large values of the e- folding time r). In this paper, we 
have presented direct observations of such young, low-mass LAEs. 

We conclude that, independently of SFR, metallicities, and extinction, the GRAPES HUDF 
LAE population is a one of very young objects with low stellar masses of a few million 
solar masses. There is no easy way to make these objects intrinsically more massive while 
keeping them intrinsically blue and without violating the IRAC flux estimates and upper 
limits shown in Table [3j We did derive somewhat larger stellar mass estimates for these 
objects when we assumed that most of their stars formed at z » 10 in our 2BP models, 
but the inferred mass of the GRAPES LAEs remained small when compared with previously 
observed Ly-a sources. 



8. Conclusion 

As part of the GRAPES slitless spectroscopic survey of the HUDF, we have identified a 
sample of faint LAE galaxies at 4.0 < z < 5.7. These objects were selected solely based on the 
detection of Ly-a emission in their low-resolution ACS slitless spectra. These spectra allowed 
us to identify LAEs down to extremely faint Ly-a luminosities (> 2.3 x 10 42 ergs s _1 ). The 
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further existence of deep space- and ground-based observations of these objects, ranging all 
the way from the rest-frame UV to the rest-frame optical, allowed us to set strong constraints 
on the masses of these objects. We found that these objects were usually best fitted by models 
with low extinction and low metallicity. While we attempted to create models using various 
SFH scenarios, both to compare our results with previous work and to attempt to maximize 
our mass estimates for these sources, we systematically derived relatively low masses. We 
conclude that the HUDF LAEs are sources with masses that are significantly lower, by several 
orders of magnitude, than previously observed LBGs at these redshifts and that they are also 
significantly less massive than previously observed LAEs at high redshifts. The GRAPES 
LAEs, with derived star formation rates of ^ 8 M Q yr" 1 and stellar mass estimates that are 
likely to be no larger than a few times 10 8 M Q , must have formed quickly, unless their SFR 
was significantly lower earlier on. The mean number density of the GRAPES HUDF LAEs, 
uncorrected for completeness, is ~ (1.25 ± 0.6) x 10~ 4 Mpc -3 . The LAE sources described 
in this paper are potentially the youngest and least massive galaxies observed to date, and 
at a time when our universe was just about as 1 Gyr old. 

This work was supported by grants GO-09793.01-A, GO-10530.11-A and AR-10299.01- 
A from the Space Telescope Science Institute, which is operated by AURA under NASA 
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tial X-ray counterparts to the HUDF LAEs using the available Chandra data, and Daniel 
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Fig. 1. — ACS G800L grism spectra of the nine Ly-a galaxies in our sample. 
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Fig. 2.— ACS, NICMOS, ISAAC and IRAC observations of GRAPES LAEs. The source is 
at the center of each postage-stamp image. ACS B 435 , V606, 1775 and z 8 so images are 1.2"- 
1.2", NICMOS Jnoo and H 1600 band images are 2.4" -2.4", and ISAAC Ks °, IRAC - 

' ilUU iWU & ' 2200 A' 3600A 

and IRAC images are 5"-5". 
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Fig. 3. — Histogram of the measured half-light radius (R50) of the GRAPES LAEs. Most 
objects are observed to have small radii, with R 50 < 1.0 kpc. 
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Fig. 4. — Colors of the GRAPES LAEs (filles squares) in z-3.6 /iin. vs. i-z space. These ob- 
jects appe ar to be extreme ly blue, similarly to previously identified high-redshift LAEs (cir- 
cles) fromlLai et al.l (20071 ). and significantly bluer than the i-dropout LBGs (open squares) 
from Eyles et al. ( 20051 ). In the left panel, colors of a z=5 single-burst galaxy are shown 
(solid line) for comparison with ages (shown in Gyr). This points to the GRAPES LAEs' 
potentially being very young, low-extinction objects. Similarly, assuming low values of ex- 
tinction (right), these sources appear to be fainter in the rest-frame optical (3.6 /xm) than 
either the LBGs or LAEs that were previously observed at these high redshifts. 
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Fig. 5. — Plots of the best double instantaneous burst model (2BP) fits (solid lines) to the 
observations (squares) and our 1 o upper limit estimates (arrows). The best, most massive 
acceptable 2BP models are also shown (dashed lines). The masses, ages, extinctions of each 
of these models are listed in Table [5j 



